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Abstract. We present results on the modelling of the ejection of a superluminal component 
in the jet of 3C 111. We propose that the component is generated by an injection of dense 
material followed by a decrease in the injection rate of bulk particles in the jet. Our model 
is supported by ID relativistic hydrodynamics and emission simulations, and is capable 
of reproducing the brightness evolution of two features, as revealed by 15 GHz VLBA 
observations. We show that other scenarios, such as an increase of the Lorentz factor in the 
material of the perturbation, fails to reproduce the observed evolution of this flare. 
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1. Introduction 

Flaring events at radio frequencies are known 
to take place in Active Galactic Nuclei (AGN), 
usually followed by the observation of new 
radio features in the parsec-scale jets (e.g., 
ISavolainen et al.l ^002). It has been shown 
that the ejection of those features, or compo- 
nents, is related to dips in the X-ray emis- 
sion fron i the active nucleus in the case of 
3C 120 (iMarscher et al.l l200 2l). and perhaps 
also in 3C 111 (lMarsched F006). The dips in 
X-rays precede the observations of new radio- 
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components. The decrease in X-ray emission 
may be caused by the loss of the inner regions 
of the disc. In this scenario, a fraction of the 
accreted material is injected in the jet and a 
new component is later observed in VLBI im- 
ages, after the material becomes detectable at 
the observing frequencies, as it evolves down- 
stream. The components are interpreted as the 
shocks produced by the ejection of denser 
and/or faster plasma in the flaring event from 
the accretion disc (Mars cher & Geaii Il985h . 
The conditions for triggering the ejection 
of the material in those radio features are 
still unknown. H ydrodynamical simulations 
dAlov et al.ll2003L A03 hereafter) have shown 
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Fig. 1. Core distance and flux density evolution 
with time of components E and F in 3C 111, 
based on the results from K08. 

that such jet perturbations produce a forward 
and a reverse structure, which would be ex- 
pected to be observed as a fast front and a 
slower back component. 

In the jet of 3C 1 1 1 (z = 0.049, 1 mas ^ 
1 pc), a very strong flaring event in early 
1996 gave rise to the ejection of two jet fea- 
tures observed at 15 GHz with the Very Long 
Baseline Array (labeled as components E and 
F -see Fig. □ and Kadleretal. 2008, K08 
hereafter). Both component trajectories can be 
back-extrapolated to similar ejection epochs 
within 3 months (around 1996.10). However, 
they show different speeds and the time evolu- 
tion of their brightness is different (see Fig.[T]|: 
whereas the inner component F is initially 
brighter (1996.82 and 1997.19) and fades out 
very rapidly (1997.66 and 1998.18), the lead- 
ing component E shows a slower decrease in 
flux density. After 1999, F has disappeared and 
E evolves accelerating. In lPerucho et"al] (l2008l) 
we propose the possibility that these compo- 
nents are the front and rear region of a single 
perturbation. Here we review this result and 
discuss on the problems that other possible sce- 
narios could face. 

2. RHD and Emission Simulations 

We have performed one-dimensional numeri- 
cal relativistic hydrodynamics (RHD) simula- 
tions in which a square perturbation in density 
is injected into a steady jet, without modifying 
the initial Lorentz factor, and relaxing the con- 
dition that the initial jet flow is reestablished 
immediately after the perturbation. We have 
substituted this by a rarefied flow, representing 
a reduction of the injection rate. In this picture, 
the original jet injection rates should be recov- 



ered after some time. However, in this work we 
only focus on the evolution of the strong ejec- 
tion and the period before the reestablishment 
of the jet flow. Multidimensional simulations 
are out of the scope of this work due to the 
computational effort required and to the one- 
dimensional character of this problem. The 
simulations have been performed using a nu- 
merical code that solves the equations of rela- 
tivistic hydrodynamics w ritten in the conserva- 
tion form, as de s cribed in iPerucho etal1(l2005h 
and lMarti et al] (Il997l) . 

The details of the simulation are given in 
the caption of Fig. |2] The top panels in Fig. |2] 
show different snapshots of the evolution of the 
square perturbation injected in a steady flow, 
in pressure, Lorentz factor and specific internal 
energy. 

Using the RHD simulations outlined above 
as input, we have computed the correspond- 
ing ID optically thin radio synchrotron emis- 
sion as seen by an observer with a line of 
sight at 19° to the jet axis (K08). For fliese 
computations, we used the numerical code and 
the procedure described in lGomez et al.l (Il997h 
and references therein. This code takes into ac- 
count all the relevant relativistic effects, includ- 
ing the light travel time delays. 

In the simulation (see Fig.|2]i, \he front re- 
gion includes the leading part of the pertur- 
bation and is identified with component E in 
K08, whereas we define the fading region as 
the rear part of the perturbation and identify it 
with component F (see Fig.|2]i. The material in 
the front region, consisting of shocked material 
from the steady jet and rarefied material from 
the perturbation separated by a contact discon- 
tinuity, shows smaller values for the pressure, 
and some acceleration due the propagation in 
the lower pression steady jet fluid. The material 
in the fading region crosses the receding rar- 
efaction that separates it from the front region 
(top panels in Fig.|2| and it is also "eroded" by 
the back rarefaction. Consequently, the front 
structure evolves increasing its size as the front 
shock incorporates material from the steady jet 
and the material from the fading region crosses 
the receding rarefaction. Thus, the front region 
consists of the forward shock structure of the 
perturbation (E in Fig.|2li, and the fading region 
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Fig. 2. Snapshots of the evolution (left to right) of a square perturbation injected in a steady jet, 
followed by a strong rarefaction. The dotted-light-blue lines stand for Lorentz factor, the solid- 
dark-blue line stands for pressure and the dashed-red lines stand for specific internal energy. The 
simulation is run with 24000 cells; the velocity of the initial flow is vy - 0.9 c; the perturbation 
is injected with a density twice that of the jet and velocity Vp - 0.9 c; the rarefied medium is 
injected after the perturbation with the same velocity as the initial flow, and pressure ten times 
smaller than that of the initial flow. Please note the change of scale in the abcissae. The bottom 
panels show the simulated total intensity emission along the jet axis at four representative epochs. 
The identification of the features in the simulation with the observed components E and F in K08 
is indicated in each panel. A jet width of 100 cells and axial simmetry is used to compute the 
emission. 



is formed by the remains of the perturbation 
that have not crossed the receding rarefaction 
(F). The synchrotron emissivity (bottom panel 
in Fig. |2]i is governed by the jet pressure and 
hence the emission evolution is very similar 
to the pressure evolution of the RHD simula- 
tions. In the emission results, the front region 
(component E) propagates without much flux 
density evolution since injection. However, the 
fading structure (component F), which initially 
shows a notably larger flux density than com- 
ponent E, rapidly decreases in emission as the 
receding and back rarefactions erode it. The re- 
verse shock (see A03) is neither relevant nor 
observationally significant in our simulations, 
as it propagates in a very rarefied medium. For 
that reason it is not shown in Fig.|2] 

Notice that the Lorentz factor values in 
Fig. 12] are those corresponding to the fluid. In 
contrast, VLBI observations provide us with 
pattern velocities. In the simulation, the veloc- 
ity of the front shock is measured to be ~ 
0.96 c (Vg** ~ 3.5c), whereas that of the fading 



region is v,. ~ 0.87 c (v" ~ 1.7c), both similar 
to those found in the observations (K08). We 
remark the fact that the velocity of the material 
in the fading region is faster than that of the re- 
ceding rarefaction (cf. Fig|2]i, as expected from 
the explanation in the previous paragraph. We 
also point out that the dilute material shown in 
Fig.|2]presents a modified velocity due to pas- 
sage through the reverse shock. 

A second simulation was performed for a 
faster perturbation, with Lorentz factor F = 
3.6, while keeping the rest of the parameters 
as in the previou s simulation. The results (see 
iPerucho et al. 2008) show that the front region 
of the perturbation is overpressured with re- 
spect to the rear region and, thus, the former is 
brighter than the latter, as shown by the emis- 
sion simulations. This is in clear contradic- 
tion with the observations of the jet in 3C 111 
(Fig. [Hand K08). The difference is due to the 
stronger front shock produced in this case. It is 
also important to mention that the wave sep- 
arating both regions is now a reverse shock, 
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instead of the receding rarefaction shown in 
Fig. |2] This is a general result for fast pertur- 
bations. 

3. Discussion 

In any of the scenarios mentioned above, 
the perturbed regions have enhanced emission 
with respect to the underlying jet. However, 
only 1) an overpressured perturbation with the 
same Lorentz factor as the underlying flow 
avoids the front region to be brighter from 
the beginning, and 2) a rarefaction behind the 
perturbation avoids the formation of a strong 
reverse shock behind the perturbation. With 
these restrictions, the second component fades 
out rapidly and then the first one dominates the 
emission, as observed for components E and F 
in 3C111. 

We note that a denser medium injected 
after the perturbation would lead to the for- 
mation of a brighter feature behind compo- 
nent F, produced by the reverse shock formed 
between the end of the perturbation and the 
medium injected afterwards. As stated above, 
this shock appears also in our simulations but 
it has a negligible effect in terms of emis- 
sion. Observational support for the inclusion of 
this tenuous material in the simulations can be 
found in the prominent emission gap following 
behind the E/F complex in 3C 111 (see K08). 
New ejection of emitting material is detected 
on the time scale of more than 2 years, corre- 
sponding to a gap width of up to 2 mas in 1999 
(cf. K08). We have exagerated the effect of the 
dilute medium in order to focus on the evolu- 
tion of the perturbation alone. However, such a 
dilute medium (a factor 10 less dense than the 
steady jet) is not needed for the conclusions of 
this work to be valid: a denser medium, but un- 
derdense with respect to the steady jet, would 
produce a similar effect, with little emission 
contributed from the reverse shock formed. 

Thus, our result explains the evolution of 
these radio-components on the basis of the sce- 
nario explained above. From this work, we can 
derive a recipe for distinguishing between ejec- 
tion events composed by denser material only 
or denser and/or faster material. On top of this, 
we want to remark that these results are based 



on the hypothesis that the injection of bulk 
particles in the jet is reduced after the injec- 
tion of the perturbation, which could be related 
to the processes taking place in the accretion- 
disk/bl ack-hole systern, in th e line of the re- 
sults in iMarscheretall (120021) . who showed the 
existence of a relation between the processes 
taki ng place in the accretion-disk and the jet 
(see iPerucho et al]|2008i for discussion on this 
issue). 

At present we are following a new radio- 
flare in this source with denser sampling. The 
aim is to perform a deeper analysis of the 
early stages of evolution of the expected radio- 
components and to compare this evolution with 
the model presented here. 
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